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We propose a method to study the second-order critical lines of classical spin-5 Ising models on 
two-dimensional lattices in a crystal or splitting field, using an exact expression for the bare mass of 
the underlying field theory. Introducing a set of anticommuting variables to represent the partition 
function, we derive an exact and compact expression for the bare mass of the model including all 
local multi-fermions interactions. By extension of the Ising and Blume-Capel models, we extract 
the free energy singularities in the low momentum limit corresponding to a vanishing bare mass. 
The loci of these singularities define the critical lines depending on the spin S, in good agreement 
with previous numerical estimations. This scheme appears to be general enough to be applied in a 
variety of classical Hamiltonians. 
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Classical or quantum spin models play a central role 
in the development of statistical physics, as they allow 
for testing fundamental principles on model systems that 
are simple enough to have a complete mathematical de- 
scription. This is particularly the case for the Ising 
model, which still serves as a toy model to develop new 
techniques both in terms of analytical methods (scries 
expansions, rcnormalization) and numerical algorithms 
(Metropolis, Swendsen- Wang) . The fact that an analytic 
solution is available in the two-dimensional case provides 
a reference point for the further understanding of criti- 
cal phenomena. While the case of the spin S = 1/2 is 
well understood, much less is known about the extension 
for spins with higher values. In the early 70's, Fox and 
Guttmann successfully developed low temperature series 
expansions, allowing for an estimation of the critical tem- 
peratures and set of critical exponents for the Ising model 
with spin S = 1 and 3/2, in two and three dimensions 
lj . This original work has later been extended to higher 
spin values 0, High temperature series expansions 
have also been used in three dimensions 0, Q and two di- 
mensions and then extended to test universality and 
hyperscaling, as well as spin-spin correlation functions 
[a |7| . Another path of investigation is provided by the 
Fortuin-Kasteleyn transformation Q , which gives a fruit- 
ful link between Ising models and percolation problems 
leading to interesting results in the continuous spin 
limit [ic| . Alternative approximation schemes have been 
proposed to treat the general spin case, such as Ising spin 
decomposition of the general spin-S model ll|, an effec- 
tive mean-field theory based on cumulant expansion fl2T | 
and Husimi tree calculation [T3|. However despite this 
activity, very few exact results have been obtained on 
the general spin-S Ising models, even in two dimensions. 



Results are available for the special case S = 1 in partic- 
ular regions of the phase diagram 14 , 15, 3, 17 1 , and for 
the case 5 = 3/2 [l8[. Recently, we proposed an alter- 
native approach to the Blumc-Capcl model where 5=1 
[l9j | . using the representation of the partition function 
with Grassmann variables [2fj| to determine the critical 
line of this model. In this article, we extend this method 
to study the general spin-S 1 model. The first step is to 
expand the partition function as a product of spin poly- 
nomials where neighboring spins are coupled. We then 
introduce a set of Grassmann variables [21( to decouple 
the spins. The price to pay for this is the lost of com- 
mutativity. However, using specific symmetries, the sum 
over the spin degrees of freedom can be performed ex- 
actly, leading to the expression of the partition function 
as a path-integral over a fermionic action on a lattice. 
Taking the thermodynamical limit of the action allows 
us to identify the bare mass of the system. Assuming 
that free energy singularities correspond to the vanishing 
mass, we obtain an excellent approximate location of the 
critical points. 

We consider now the general Hamiltonian on a two- 
dimensional (2D) lattice of size L x L 
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where are Ising coupling constants and Ao a split- 
ting crystal field that favors small spin values. It can 
also represent a chemical potential in the Blume-Capel 
model. In the Ising case, u^ nn = 1, and therefore Aq 



2 



does not play any role. This crystal field term can be re- 
placed by any potential V(a mn ) depending on the square 
of the local spin. The spins a mn can take 25+1 values 
with a mn = — 5, — 5 + 1, • • • ,5. We assume here for sim- 
plicity that Ji = 3% = 1 and K = (iJ\.i = 1/fcsT the 
inverse of temperature. The partition function, which 
represents the sum over all possible spin configurations 
Z = 1t{ CT } exp(— 0H), contains products of the Boltz- 
mann weights exp(Kaa') (where a and a' are neighbor- 
ing spins) which take q + 1 = 5(5+ 1) + 1 different values 
if 5 is an integer, and q + 1 = (5 + l/2)(5 + 3/2) values 
if 5 is half-integer. Since there are q + 1 possible values 
for each Boltzmann weight exp(Aercr'), we can project 
each of them onto a polynomial function of degree q in 
variable aa': 

q i 
exp(Aercr') = ^ u k (aa') k = u (1 + x a aa'), (2) 

fe=0 a=l 

where the 5 + 1 constants are determined by solving 
the linear system of q + 1 equations satisfied by the 
above relation. Note that the demonstration below holds 
for any function of aa' in the Boltzmann factor, in the 
case where the Hamiltonian (JT]) includes quadrupolar 
interactions {aa') 2 for example fl8| . Only the coeffi- 
cients Mfc are different. In the Ising case, 5=1/2 and 
q = 1, we have exp(Acrer') = ch(A:/4) + 4sh(A74)crcr', 
and it = ch(A"/4), ui = Ash(K/A). In the Blume- 
Capel model q = 2 and it is easy to show that 
uq = 1, m = sh(A) and «2 = ch(A") — 1. For 5 
integer we always have uq = 1, and from Eq. 

Ul<fe< 9 = MQ E ai<ai <-<« t EaiZaa • ■ • a 



We also set 



Uk>q+1 



0. 



Let us now introduce q pairs of Grassmann variables 
{aZ\ni®mn) a * eacn s i te f° r * ne horizontal direction and 
q additional pairs (b mn ,b mn ) for the vertical direction 



SHU. Here a takes the values 1, . . . q. In total, there 
are Aq Grassmann variables on each site of the lattice. 
The Ising model is usually represented by two pairs of 
Grassmann variables per site which can be reduced after- 
wards to one pair [2(| ■ We then use the following integral 
representation for each couple of terms (1-\-x a a mn a m +i n ) 
and (l+Xaamnamn+i) that appear in Boltzmann weights 
Eq. ©: 

1 + ^a^mn^m+lrt 

dd mn da mn e mn (1 a mn IJ, mn){\ + X a d, mn a m -^i n ) , 
1 + Xa&mn&mn+l = (3) 
^mti^mti 6 mn m ™(l + b mn a mn ){\ + X Q 6 mn cr mn -|_i). 



partition function can be written as 
Z 
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Tr 
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where A = — (3Aq. Notice that inside the integral sym- 
bol, the pairs of link factors in brackets (A™ n A° +ln ) 
and (Bmn^mn+i) can move freely. In particular, we can 
rearrange the products over a and put together link fac- 
tors having the same site indices (m, n) using the mirror 
ordering symmetry [20| 




m+ln I ' 



where the arrows indicate that the product is ordered, 
i.e. increasing label a in the first product from left to the 
right and in the second one from right to the left. For con- 
venience, WC SCt O mn = U. a A mm O m+ln = Ua^m+ln 
and V mn = ria-Smn, Ann+1 = TllCn the 

partition function can be rewritten as 
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At this stage, we use the mirror and associative symme- 
tries which were applied to the Ising model [22j and which 
are still valid here to rearrange the operators O, 0,P and 
■p. The computations are until now identical to the Ising 



case treated in references 22, 231, in the sense that we 



obtain an expression of the partition function with a set 
of anticommuting operators we would like to rearrange 
in order to perform the sum over the individual spins. 
The only difference is that the previous operators O, O, 
V and V are in general more complicate functions of the 
Aq Grassmann variables coming from the decomposition 
given by relation ((2|) . In principle boundary terms should 
be treated separately in order to obtain the exact finite 
size partition function. Periodic boundary conditions can 
be treated rigorously for finite lattice 0,[i3,[23 but this 
is inessential in the thermodynamical limit L — > oo we 
are interested in here. We consider instead free boundary 
conditions, leading to the exact expression : 



From these expressions, we introduce the link factors 
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Under this form, the spins can individually be summed 
over from <tl„ to o\ n for any given n. We introduce the 
following weights W mn which include all the dependence 
on the individual spin a mn 



<r ra „=±l 
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where we have defined the following sets of Grassmann 
variables 

x 9 _ 0+ io^ l _ 1 „ if a = 1, • • ■ ,q, 
x 2q - a +ibl q l ~"i 1 if a = q + 1, • • • , 2q, 



a— 2q 
^ran 



if a = 2g+ 1, 
if a = 3q + 1, 



,3g, 
,4g. 



The partial sum {5) can be performed by noticing 
that only products involving an even number of o~ mn 
give a non-zero contribution. We also define a k 
J^ = _ s cr 2fe exp(Acr 2 ) and the ordered products q 



(fe) _ 



E« 1 <»2<-<« i CC • • • Cn with 9mn = 1- With these 
notations, it is easy to show that the partial sum ([5|) gives 
the commuting objects 



2./ 



Wrr 



ran ' 



fc=0 



with qii n ' = c mn ■ ■ ■ c^ in the term of highest de- 
gree in Grassmann variables. Finally the fermionic 
representation of the partition function reads Z = 
ul L ~ %x {aab ^Y[ mn W mn . This is the exact fermionic 
representation of the partition function for any given 
spin-S' model. 

In some cases, the weights W mn can be easily exponen- 
tiated. For the Ising model, the argument of the expo- 
nential is purely quadratic in the c mn 's and therefore the 
partition function can be written as a determinant [23j |. 
In the Blume-Capel model, the argument of the exponen- 
tial is a polynomial of degree 8 in Grassmann variables 
since there are 8 independent variables (Aq — 8) In 
general we expect the argument to be at most a polyno- 
mial of degree 4g in these variables, which can sometimes 
be reduced by partial integrations. Except for the case 
q = 1 the partition function is not solvable. 
In the thermodynamical limit L — > oo, however, we ex- 
pect to be able to identify from the effective theory a 
massive and pure kinetic contributions in the infrared 
region where the continuous momenta k are small. The 
condition of criticality is determined usually by the van- 
ishing mass m of the effective theory. For example, in 
the Dirac or Majorana representation of the Ising model, 
the free energy, which is the integral over the Brillouin 
zone of momentum-dependent quantities In (771 2 + k 2 ), is 



singular at m — 0. The action is determined by the ex- 
ponentiation of the W mn quantities, depending on the Aq 
Grassmannian fields, and is made of a local part, con- 
taining all the local interactions, including 2-fcrmion, 3- 
fermion etc.. interactions at a given site, and a kinetic 
part, containing all the terms involving space deriva- 
tives of different orders. It is difficult to obtain the 
full fermionic action with all the kinetic terms in the 
general case. In this paper, we will neglect the latter, 
assuming that their contribution by renormalization to 
the mass is negligible near the critical point. This is 
true for the Ising model where the kinetic part is purely 
quadratic in Grassmann variables and does not renor- 
malize the mass. Space symmetries of these derivative 
terms also could prevent any renormalization. Then we 
will show that the contribution to the partition function 
coming from the local part only can be computed ex- 
actly, which is not a quadratic action but a polynomial 
of degree 4<j, the number of variables involved, and this 
defines a bare mass. We need for this to define first the 
formal derivatives of Grassmann variables (26l , for exam- 



ple: d x a n 



.in and d v a r , 



Then the c's coefficients can be expressed in term of 
these derivatives such as c mn — x q (a^ nn — dxd^). In 
the limit of large L and in the Fourier space, the first 
order derivatives account in the action for a small con- 
tribution in momenta k = 27r(m, n)/L, with m, n <C L 
positive integers, when amplitudes |k| become small. In 
this infrared regime, we assume here that we can neglect 
the derivatives : c mn ~ XgO^, • • • , c mn ~ xxa x mn and 

c rnn - ^mn, '",4 - X l b lan ■ the Weights W m n are 

then all decoupled, and the following bare mass ms can 
be defined: 



m s 



[ ]J da mn da a m ndb a mn db a mn e a 



W mn .(6) 



The different integrals in Eq. © can be evaluated ex- 
actly by noticing for example that the arguments of the 
exponential bf ml bf nn can be combined with a a mn {x a d mn ) 
that appears in some of the q^ products to give a con- 
tribution x a . Indeed using the Grassmann integration 
rules J da. a = 1 and J da.l = 0, we can write 



da a da a db a db a e a " 

UjUl mn UjU 'mn UlU mn uu van c - 



i( x a& m n) x a ■ 



Since the q( 2k ^ are ordered, there are also signs to take 
into account and coming from moving the variables c mn 
before integration. We obtain after some combinatorial 
algebra 
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2q 
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(7) 
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TABLE I: Critical temperatures at Ao = 



Spin S 


5 = 1/2 


5 = 1 


5 = 3/2 


5 = 2 


g 


1 


2 


5 


6 


t c 


2.269 185 


1.673 971 


1.456 694 


1.337 812 


Refs. 


2.269 [6] 


1.689 [1], 1.695 [27] 1.461 [6, 28, 29J 


1.336 [6] 






1.694 [6], 1.681 [28] 






Spin S 


5 = 5/2 


5 = 3 


5 -> oo 




9 


11 


12 


00 




tc 


1.262 542 


1.210 534 


0.925 148 




Refs. 


1.257 [6] 


1.203 [6J 


0.915 [6, .10] 





where we have defined the following quantities Rq = Uq, 

k 



Rk = ^2 UiUk-icr(l, k - I), 



(8) 



i=o 



and u(k, I) = 1 if k and / are both even, and cr(fc, I) = — 1 
otherwise. We now apply this result to different cases. 
For the Ising model (S = 1/2) we obtain mi/2 = 
2e A / 4 [l — sh(if/2)], which vanishes at the normalized 
Ising critical temperature t c = T c /S 2 = 2.269 185, in- 
dependent of Ao. For the Blume-Capel model (S = 1) 
we find mi = 1 + 2e A [l - sh(2K)} and for S = 3/2: 

m 3/2 = 2e A / 4 [l - sh(K/2)] + 2e 9A / 4 [l - sh(9X/2)]. (9) 



For general spin 5, we can show that 



m s = J2 e Aa2 [l-sh(2a 2 A')]. 



(10) 



a=-S 



Equation (fTTJj) gives the expression of the bare mass of 
a general spin-5 system, taking into account all possi- 
ble local fermion-fermion interactions. To go further we 
propose to extract some physical information from the 
previous result in different cases. In particular, for the 
simplest non-integrable case (Blume-Capel model S = 1), 
it is also possible to write explicitly the fcrmionic action 
and to check Eq. (JTUJ) fill]. Note however that even in 
this case a vanishing mass is a necessary but not suffi- 
cient condition to have a critical point : kinetic terms of 
higher order that appear in the effective fermionic action 
can change the nature of the singularity. For the Blume- 
Capel model for instance the critical line terminates at a 
tricritical point which can not be predicted by the mass 
alone 0. 

Tabulated values of t c at Ao = are given in Table |T] for 
several S, and compared with numerical results (Monte- 
Carlo simulations, high-temperature expansions) given in 
the literature. In general, the agreement is good. For in- 
teger values of S (Blume-Capel model), the critical line 




t=T/S 

FIG. 1: Critical line for different values of 5 = 1,2,3, 
5 = 3/2, 5/2 and 5 = oo. In the latter case, the critical tem- 
perature at Ao = is t c — 0.925 148, and at t = the curve 
reaches the solution Ao = 4/\/3. The inset is a zoom on the 
region Ao > 2. The critical curve for 5 = 3/2 has an asymp- 
tote (dotted vertical line) at t = 2/91og(l + \/2) ~ 0.252 132. 



goes from the Ising critical value t c = 2/ log(l + v / 2) when 
Ao — * — oo to the terminating point (t c = 0, Ao = 2) 
continuously. For half-integer values of 5, there exists 
in general an asymptote in the (t = T/S 2 , Ao) plane. 
Indeed, for S = 3/2, Eq. ^ predicts the solution 



Ao 



9/ 



log 



1 - sh(2/9t) 



1 - sh(2/i) 



(11) 



which is bounded by t c = 2/91og(l + y/2) ~ 0.252131 
below which there is no second-order critical line (see 
FigO}. In the large integer S limit, the model defined 
in Eq. |T]) is described by a continuous variable — 1 < 
Xmn = Pmn/S < 1 (continuous Ising model). We can 
obtain the limiting value of the mass, Eq. (|10p becoming: 



ms>i ~ 5v 2t 



li It 



dx e - Aox2/2 



l-sh(a; 2 ) . (12) 



We observe that the rescaled mass mg/S vanishes 
when A = at t c ~ 0.925 148, in fairly agreement 
with numerical works @, [Toj ] . and there is a non 
trivial solution at t = which is simply given by 
A = 4/\/3 ~ 2.309 401. Contrary to the finite S 
models, the critical field predicted here is not equal to 
Ao = 2, which gives the location of the first-order transi- 
tion at zero temperature, but takes a slightly larger value. 

By using Grassmann algebra to represent the parti- 
tion function of spin-5 Ising Hamiltonians on 2D square 
lattices as a fermionic theory, we were able to obtain 
the exact expression for the bare mass of the action, in- 
cluding all the possible local fermionic interactions. This 
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result gives at least precise though approximate location 
of second-order critical points in the (T, Ao) plane. This 
scheme and its main consequences, formulas (J7]) and ((HJ), 
are general enough to be applied in a variety of classi- 
cal Hamiltonians with next-nearest neighbor interactions 
and crystal field like potentials, and possibly Potts-like 
models as well, with a suitable choice of polynomial rep- 
resentation of the Boltzmann weights. 
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